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Metallic thin film ferromagnets generically possess spiral states that carry dissipationless spin 
currents. We relate the critical values of these supercurrents to micromagnetic material parame- 
ters, identify the circumstances under which the supercurrents will be most robust, and propose 
experiments which could reveal this new collective transport behavior. 
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In ferromagnetic metals and semiconductors quasipar- 
ticle states can be manipulated by external magnetic 
fields that couple to the spin-magnetization-density col- 
lective coordinate. This property is responsible for re- 
lated robust magnetoresistance efi'ects that occur in var- 
ious geometries such as anisotropic magnetoresistance in 
bulk samples, giant magnetoresistance in metallic 
multilayers, and tunnel magnetoresistance ||^,^ in tun- 
nel junctions. In this paper we propose a distinctly dif- 
ferent type of spin-dependent transport effect, in which 
spin current is carried collectively rather than by quasi- 
particles. Because the spin current is non-zero when its 
quasiparticles are in equilibrium, it is carried without 
dissipation. This spin-supercurrent state occurs only in 
easy-plane ferromagnets and will be robust only when 
anisotropy within the easy plane is weak. We propose an 
experiment to observe this effect in thin films of ferro- 
magnetic metals. 

The key observation that motivates this proposal arises 
by considering the class of excited states obtained from 
the ferromagnetic ground state by following its adiabatic 
evolution as equal and opposite constant vector poten- 
tials are introduced for up and down spins, with the spin- 
quantization axis perpendicular to the ferromagnet's easy 
plane. The many-particle Hamiltonian is 
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with 7Yci-oi — J2i<j ^'^ ^ The vector potentials 
for spin a are A^. = Qa{^c/e) with Q| — — Qx = Q/2. 
They can be removed by gauge transformations, multi- 
plying the many-particle wavefunction by exp(zQ(j-r) for 
each electron with spin cr. In a paramagnet, the ground- 
state wavefunction would, therefore, evolve trivially with 
Q and the ground-state energy would be independent of 
Q. In a ferromagnet, however, a change in the phase 
relationship between up spins and down spins alters the 
magnetic order and will change the energy. 

We start with a ground state that has a spontaneous 
spin density along the x direction. Its magnitude is 



m{r) — (^'|(r)^'|(r))o, where ^^(r) is an electron field 
operator for an electron with spin a |, and (...)o de- 
notes a ground-state expectation value. For small Q, 
^'cr(r) — > exp(iQcr • r)\E'cr(r) and the order parameter ro- 
tates in the x — y plane as a function of r, 



(s(r))Q = mQ[cos(Q • r)x + sin(Q • r)y] 



(2) 



i.e., it forms a spiral spin state. As Q increases, the order 
parameter's spatial dependence will ||] cause a decrease 
in the magnetic condensation energy and in the magni- 
tude of the order parameter. Dependence of the ground- 
state energy E on Q implies that these many-particle 
eigenstates have finite current densities with equal mag- 
nitude and opposite direction for up and down spins, 
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where e(Q) is the energy per unit volume. 

As we discuss at greater length below, easy-plane 
anisotropy ascribes a topological character to the 
wavevector Q, so that these currents can decay only by 
phase slip processes that have large barriers, i.e., these 
are dissipationless supercurrents. 

Equation (^) is similar to the connection between the 
exchange-coupling of ferromagnets separated by a tunnel 
junction and the spin currents that flow between them 
Our proposal for supercurrents in ferromagnets is related 
to Anderson's discussion of superconductivity Q in terms 
of magnetic order in an effective spin model; the physics 
of the two ordered states appears similar if a particle-hole 
transformation is made in one of the spin subspaces. The 
supercurrents we propose are also related to those sup- 
ported in double-layer quantum Hall systems l8|-[T2|?], 
where ordered states form that are describeable either as 
pseudospin ferromagnets or as electron-hole pair conden- 
sates [ p"3|jr^ ]. In fact, the role of magnetic anisotropy in 
controlling the observability of these supercurrents is con- 
nected in part with the role of band hybridization terms 
in controlling the observability of collective electron- hole- 
pair transport in excitonic insulators p"^-p^. 
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To illustrate these ideas we consider the simplest pos- 
sible microscopic model of a metallic ferromagnet, a 
fermion gas with delta-function repulsive particle-particle 



interactions U5{ri — rj) treated in a mean- field approxi- 
mation. The unrestricted Hartree-Fock Hamiltonian for 
the ordered state with wavevector Q = Qx is 



where ek = h'^k^/{2m), the ordered-state quasi- 
particle energies are E]^^± = [ek+Q/2 + ek-Q/2 ± 
V(ek+Q/2 - + 4/i2]/2, and the effective mag- 

netic field which splits the quasiparticle bands by ±/i 
with h ~ U niQ is determined sclf-consistently by 



U V(ek+Q/2-ek-Q/2)^+4/i^' 

The prime on the sum in Eq. indicates restriction to 
those wavevectors for which only the lower-energy quasi- 
particle state is occupied. 

The procedure described above provides a mean- 
field approximation to the ferromagnet spin-supercurrent 
states. In Fig. we plot quasiparticle bands for a typ- 
ical model of this type at a moderately large value of 
Q = Q.^kp, where kp is the Fermi wavevector for zero 
order parameter, and ef is the Fermi energy |jT^. The 
product UD^ep) was chosen to be close to experimental 
values for Co and Fe (taken from Ref. [0). In this case 
the Stoner criterion |2l| for mean-field ferromagnetism, 
UD{ep) > 1, is satisfied. In Fig. ^ we plot the order 
parameter toq, the magnetic condensation energy econd, 
and the spin supercurrent density j as a function of the 
ordering wavevector Q. Note that the current density is 
proportional to the derivative of the condensation energy 
in agreement with the more general discussion above. 

Our calculations demonstrate that spin supercurrents 
are possible in states with equilibrium quasiparticle pop- 
ulations; elastic scattering from occupied to unoccupied 
states cannot provide the current decay mechanism famil- 
iar from the standard theory of metallic transport. To 
establish the stability of the spin currents it is, however, 
still necessary to show that the spin-supercurrent state 
is stable against infinitesimal distortions of its order- 
parameter field. In what follows we demonstrate that 
magnetic anisotropy is necessary for stability. Since real 
metallic ferromagnets are much more complex than the 
toy model system discussed above, we now turn to a phe- 
nomenological approach that will allow us to relate crit- 
ical currents to known micromagnetic parameters. 

We consider a generalized Landau-Ginzburg model 
for the dependence of an easy-plane ferromagnet 's free- 
energy density |po| , p^ on its magnetic state: 

/ = -|a|M • M + ^(M • M)2 + i|VMp + KAll (6) 

The free energy of this model is minimized by a constant 
magnetization in the x — y easy-plane with magnitude 



Mq — \/\a\/p. The resulting dependence of energy den- 
sity on magnetization orientation at fixed magnitude al- 
lows us to identify A = AAIq with the exchange constant 
and K — KMq with the uniaxial anisotropy coefficient of 
the ferromagnet 's micromagnetic energy functional [ pO[ . 
The spin-supercurrent state has magnetization MqIx) — 
MQ{cos{Qx),sm{Qx),0), where = {\a\ - AQ'^)/(3 is 
decreasing with Q as in our microscopic calculations. Us- 
ing Eq. (^), we find that the spin supercurrent density is 

j(Q) = ^(H-iQ^), (7) 

reaching a maximum at Qph where AQ^f^ — \a\/i. 

Expanding around the spin-supercurrent state free en- 
ergy extremum, we find that 

5,f = 2(iMlMl + A\VMa\^ + i|VMp,,|2 

+2QA{Mad,Mph - Mphd^^Ma) 

^{k - AQ^)Ml^ A\S/M,\'^, (8) 

where Ma and Mpu are the amplitude and phase fluctu- 
ations of the easy-plane magnetization (the projections 
along and perpendicular to Mq), while Mz is the hard- 
axis fluctuation. The translationally invariant kernel of 
this quadratic form has three wavevector (p) dependent 
eigenvalues: 

K± = /3M| + ± ^j3^M^ + AA^Q-^pl (9) 
Kzz^k-AQ^ + Ap'. (10) 

It follows from Eq. ^ that the spin-supercurrent state 
is stable against easy-plane fluctuations provided that Q 
is smaller than Qph', at larger values of Q, energy can 
be lowered by phase separation into regions with larger 
and smaller Q. For the soft ferromagnets we have in 
mind, however, it is the out-of-plane fluctuations, de- 
scribed by E q. (p^ , that become unstable flrst. For 
Q > Q, = y^K/A, the spin supercurrent can relax by 
tilting out of the easy-plane to one of the poles and un- 
winding phase with no energy cost. In Table ^ we list 
Qz values and the corresponding critical current densi- 
ties jcrit = j{Qz) for some common soft thin fllm mag- 
nets, including only the shape (magnetostatic) contribu- 
tion iiTshapc = Mo-^o /2 to K. From our model calculation 
and the results shown in Fig. p| we conclude that Qph is 
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typically of the order of the Fermi wavevector kp, i.e., 
much larger than Qz- To estimate jcrit we can, therefore, 
linearize Eq. (|^) in Q, 

jcrit = 2|^Q^ = 2^y/AKsh^p^ , (11) 

to obtain the large critical currents listed in Table |[ 

We have so far neglected magnetocrystalline 
anisotropy, since it is much weaker than shape anisotropy 
in the situation we have in mind. It does, however, break 
rotational symmetry within the easy plane and has the 
tendency to fix the phase and, thus, to suppress the 
supercurrents. When an in-plane anisotropy term is 
included in the energy-density functional, extrema at 
small phase winding rates consist of weakly coupled soli- 
tons in which the magnetization goes from one in-plane 
minima to another. (Q = ONg/L where 9 is the angle 
between in-plane minima and Ns/L is the soliton den- 
sity.) The energy density at small Q is proportional to 
the number of solitons. As a consequence, the minimum 
spin-current density jmim that can be supported by a 
spin-supercurrent state is non-zero. To estimate jmin for 
cubic materials we include the leading-order bulk cubic 
anisotropy in the energy density, k['^^ sin^ ip cos^ tp 
where ip is the angle of the order parameter within the 
easy plane. For small Q the functional is minimized 
by a kink soliton solution. By evaluating the energy of 
in-plane solitons of this model, we find from Eq. (||) that 

ii^^^Jj^ (12) 

jcrit 47r y Xshapc 

which is of the order of 1.5% (see Table |). (100) hep 
Cobalt thin films with in-plane easy-axis will typically 
have still smaller values of jmin because of the higher 
hexagonal symmetry. From these considerations, we con- 
clude that spin supercurrents will be observable at mod- 
erate current densities only in materials that have weak 
magnetic anisotropy within the easy plane. Because 
of their extremely weak magnetocrystalline anisotropics, 
homogeneous permalloy samples might be ideal candi- 
dates for the experiments proposed below. Although the 
course grained in-plane magnetic ansiotropy can in prin- 
ciple be fine-tuned to zero, spin-rotational invariance in 
the easy-plane will always be broken by disorder terms in 
the microscopic Hamiltonian. Since dissipationless spin 
supercurrents will not occur if these disorder terms are 
too strong, the effects we propose are more likely to be 
observable in homogeneous alloys. 

One possible experimental arrangement in which this 
collective transport phenomena could be detected is illus- 
trated schematically in Fig. |^. An easy-plane thin film 
ferromagnet (F) is connected to four spin-selective leads 
(full spin polarization in the leads is optimal but not 
required) that feed opposing up and down spin currents, 



where "up" and "down" refers to the direction perpendic- 
ular to the thin film. We emphasize that even with recent 
advances in transition metal ferromagnet spintronics, re- 
alizing a system with this geometry represents an exper- 
imental challenge. In this setup, a quasiparticle current 
would flow dissipatively between upper and lower leads 
on both the left and right hand side of the thin film ferro- 
magnet. A sizeable voltage drop (measured, e.g., between 
the upper leads), proportional to the injected currents, 
would result. Its exact value depends on the resistivity 
of the ferromagnet and on details of the geometry and 
is not a concern here. If the collective transport effect 
predicted here occurred, however, currents with oppo- 
site spin would flow without a voltage drop across the 
sample, from left to right and vice versa. Dissipationless 
current flow in the bulk could still be masked by resis- 
tance in the film-lead contacts or by collective spiral wave 
phase-slip processes. Our uncertainty in the magnitude 
of the contact resistances compared to the quasiparticle 
resistance makes our proposal somewhat speculative. A 
collective element to the spin transport could be unam- 
biguously identified by driving the critical current density 
j through either the maximum or the minimum current, 
jcrit or jmin, Or by reversing the spin orientations of the 
leads on one side of the sample. The later change would 
have no effect on the measured voltage if the current were 
carried entirely by quasiparticles but would increase the 
voltage if part of the current was carried collectively. 

In conclusion, we have examined circumstances under 
which dissipationless spin supercurrents, associated with 
spiral magnetic order, can occur in thin film ferromag- 
nets. We have estimated critical values of these supercur- 
rents and proposed an experiment to generate and detect 
this new collective transport behavior. 
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FIG. 2. The order parameter mg normalized to elec- 
tron density rie, the magnetic condensation-energy density 
Econd normalized to the energy density of the disordered 
state, eo = (3/5)neeF, and the spin supercurrent density 
j = J I — — j| normalized to jo — enehkp /m as a function 
of the ordering wavevector Q for UD{eF) = 1.5. The dashed 
lines indicate an instability regime against phase separation 
into regions with larger and smaller Q. 
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FIG. 1. Quasiparticle bands Ek,+ (upper solid curve) and 
Ek.- (lower solid curve) for Q = 0.5fcF, UD{eF) ~ 1.5 and 
ky = kz = 0. For comparison we also show the dispersion 
ek+q/2 and ek-q/2 for zero order parameter (dashed lines). 
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FIG. 3. Schematic illustration of one possible experimental 
set up to prepare a spin-supercurrent state. 





Fe 


Go 


Ni 


MoMo [T] 
4[pJm-i] 
Ki"^ [MJ m-3] 
Q. [nm-'] 
jcrit [A cm"^] 

Jmin/ J?crit 


2.15 
8.3 
0.048 
0.47 
1.19 X 10^ 
0.013 


1.81 
10.3 

0.36 
1.11 X 10** 


0.62 
3.4 
-0.005 
0.21 
2.2 X 10* 
0.015 



TABLE I. Saturation moment /loMo, exchange constant A, 
cubic anisotropy constant K['^\ critical wavevector Qz, criti- 
cal spin current density jcrit , and ratio of minimum to critical 
spin current density for common soft thin film ma gnet s. The 
values for fioMo, A, and k['^^ are taken from Ref. [pO[. 
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